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I Abstract: We study the properties of the electroweak sphaleron on a finite 

lattice. The cooling algorithm for saddle points is used to obtain the static 
O ■ classical solutions of the SU(2)-Higgs field theory. Results are presented for 

Mh = oo, Mw, jMw- After performing finite size scahng we find good agree- 



X 



ment with the results obtained from variational approaches. Of relevance for 
numerical determinations of the transition rate is that the lattice artefacts 
are surprisingly small for Mw ~ Mh- 



> 

O 1 Introduction 
(N 

ly-^ ! In this paper we will study the sphaleron solutions for the SU(2)-Higgs field theory, using 
0^ I the lattice approximation and an algorithm to find saddle-point solutions. The sphaleron is 



a solution of the static equations of motion, i.e. a stationary point of the energy functional, 
which has precisely one unstable direction. This direction corresponds to the tunnelling 
path associated to the (approximate) instanton. Due to the spherical symmetry, variational 
r-| ! analysis using a radial ansatz has provided accurate results quite some time ago [|l], ||. 

However, due to the recent interest of studying the sphaleron transition rates on a lattice 0, 
the question arises how big the lattice artefacts are for the particular sizes of lattices that 
are employed in the numerical analysis. The lattice destroys the rotational invariance and 
a variational analysis does no longer seem very practical. Furthermore, in the absence of 
rotational symmetry in the continuum, the method discussed can be used with the same 
ease. 

We have reported earlier on the sphaleron solutions where the length of the Higgs 
field is frozen. In the unitary gauge this means that we only need to consider gauge degrees 
of freedom. We recall that above Mh = 12Mw the sphaleron undergoes a series of bifurca- 
tions acquiring at each bifurcation an additional negative mode, while new solutions, 
so-called deformed sphalerons split off. For infinite Mh, where the model is identical to the 
gauged non-linear sigma model, there is an infinite number of solutions ranging in energy 
from bAlMw/aw to the energy of the lowest deformed sphaleron b.OTMw/aw, which has 
only one negative mode (the number of unstable modes increases with increasing energy). 
These solutions are related to the electroweak skyrmions H. 
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Here we will include the scalar field in the analysis to allow study of the electroweak 
sphalerons (at 6^ = 0) for a more interesting range of parameters. We will report results 
for Mff = M\Y and = jM^, the latter value corresponding to Mh ~ 60GeV, the 
present experimental bound for the Higgs mass 0. Since for finite values of the Higgs 
self-coupling the scalar field is allowed to vanish at the center, these solutions are smoother 
(have smaller lattice artefacts) than for the electroweak skyrmions. We first present the new 
algorithm to find the extremum of the energy functional, based on minimizing the square 
of the equations of motion. A careful analysis of the finite size scaling is performed, to 
allow for a reliable extrapolation to the infinite volume limit. The results agree accurately 
with those obtained from the variational analysis. For = My/ the lattice artefacts are 
to a good degree described by the formula E = Eq — 0.3(aMvi/)^ — 0.3{aM]yy, whereas 
the volume corrections are described by 3.641 + 18.1{MwL)~^e~^'^^^ (the infinite volume 
variational result p is 3.6417) all in units of Mw/ctw, where aw = 5'^/47r is the electroweak 
fine-structure constant. 

2 The model 

The dynamical variables for the SU(2)-Higgs model on the lattice are the gauge group 
variables V^(a;), defined on the link that runs from x to x + fi, and the Higgs field in the 
fundamental representation of SU(2) (a complex two-component spinor) defined on the site 
X. This Higgs field can be represented by its length p{x) (in the continuum this neutral 
Higgs field will be denoted by (j){x)) and a SU(2) matrix a{x), which is associated to the 
gauge degree of freedom and can be reabsorbed into the links via the change of variables |§] 

Uf,{x) = a{x)Vf,{x)a{x + fi). (1) 

This gives the Higgs model in the unitary gauge. The lattice action is (?7^(a;) = ^y;) 

9 X { /1,U X fl' fj. ) 

(2) 

For n = A [n = 3) the continuum action (energy) functional is recovered by rescaling the 
fields and coupling constants. Introducing a lattice spacing a, to convert to dimensionful 
parameters, one first scales the fields to get the correct normalizations for the kinetic terms. 

2 2 

U,{x)=exp{aA,ix)), A,{x) = -tgA^ix)^, p\x) = ^<p\x), (3) 
where are the Pauli matrices. The continuum parameters Myy^ Mh and A are given by 

A = |^, {aMwf = —, {aMHY = ^, (4) 
with V the lattice vacuum expectation value 



Introducing the parameters 

R = 2M^a^, r = Mh/Mw, (6) 
one can eliminate A and k, in favour of these more physical parameters 



r^K^ -(32 - Rr^) + J{32 - Rr^f + 16r^ 
^ = ^' = 2? • 

Note that for r oo, f — 1 and n ^ R. 

In this paper we are interested in the energy functional, with Uo{x) = I, and all fields 
time independent. Note that restricting the sums over the indices to three dimensions 
leaves an extra term —2np'^{x) from the time component of the hopping term. We have 
chosen our conventions such that the gauge coupling constant can be factored out, allowing 
us to express the energies in units of Mw/oiw 

8 = E{ETr(l-"n)-KEp(a:)p(a:+/i)Tr([/,(x)) + (l-2K)/(a;) 



+X{p'{x)-iy-Co}. (8) 



From now on all indices are assumed to run over the values 1-3. The constant Cq normalizes 
the vacuum (f/^(a;) = / and p(x) = v) energy to zero, 

Co = {l-8Ky + X{v^-lY = \-\{aMH)\aMwf. (9) 

For ease of reference we quote the continuum expression for the energy functional in the 
unitary gauge using our conventions {F^y = d^Aj,-duA^ + [A^, A^] = -igF^j^Ta/2) 



^ = ^1 d'x\Tr {F^,y\+ I (fx 

3 Cooling 



(9,0)^-1 Tr (4)02 + A (ct?-Ml/st' 



(10) 



Cooling algorithms are designed to find a solution for the equations of motion associated 
to a local minimum of the energy functional. It is relatively easy to write down the lattice 
equations of motion. In particular it should be noted that the energy functional depends 
linearly on the links. One finds 

du,ix)S oc U,{x)Ul{x)-U^{x)Ul{x) = 0, d,^x)S oc p{x){l-2K+2\{p\x)-l)]-p{x) = 0, 

(11) 

where ^ _ _ M x 

u,{x) = iKp{x)p{x+p) + u,{x-Q), ^/.(^;0) = E(-n + "U)' ^^2) 

p{x) = \kY,p{x+p)1:i{U^{x))+ p{x-p)i:i{U^{x-p)). (13) 

The equations of motion for the links are solved by 

U,{x) = ±iJ,{x)/\\iJ,{x)\\. (14) 
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The positive sign is to be taken in order for the solution to have a smooth continuum 
hmit. The solution for the scalar field p{x) is given by the root of a cubic polynomial. 
If 1 — 2/? — 2A > 0, equivalent to the condition {aMnY < 12, there is only one real root, 
p{x) = ps{p{x)), where 



Ps{p) = 



1-2k-2X\ ^ p 



?.2 



6A 



16A2 



+ 



4A \ 



1-2k-2A\ ^ p 



n2 



6A 



16A2 



(15) 



Cooling is performed by iterating these equations, i.e. replacing the link and the scalar 
field by the right-hand side of these equations, sweeping in a particular order through 
the lattice. With only nearest-neighbour interactions, checkerboard-type updates are most 
efficient and allow for vectorization of the algorithm. We use this cooling to first bring a 
random configuration down to one that is smooth. But since the solutions we are interested 
in have an unstable direction, we should switch to an algorithm that does not make the 
solution decay along the unstable direction (to the vacuum). This is achieved by taking the 
square of the equations of motion as the minimizing functional [|lOl, and devising an efficient 
algorithm for minimization []TT], ^. There are of course more sophisticated algorithms to 
avoid decay along an unstable direction, but they tend to require information on the Hessian 
of the energy functional, which is expensive for large lattices. 



4 Saddle-point cooling 



We define S by summing the squares of the equations of motion, {du^{x)£Y ^'^^ {dp{x)£Y-, 

^=^E{Tr(^MW^JW-[f^MW^JWF)+^('^^^) [l-2«:+2A(p2(x)-l)]-p(x))'}, 

(16) 

where / is an arbitrary positive constant. One can show that in the continuum limit 



S 



9>-4A0=^ + iTr(A^)0+iMi0 , (17) 



which has the dimension of M^r. Consequently, we will quote values of S in units of 
M^/aw- For r finite, k has a non-zero limit when a — (e.g. n{a = 0,r = 1) = 0.1245), 
we therefore took / = 1. Saddle-point cooling introduced in ref. |jlT[ is designed to minimize 
5* down to its minimal value of zero. The value of S* is a direct measure for how close the 
cooled configuration is to an exact lattice solution. 

Finding an algorithm to minimize S is more complicated due to the quadratic depen- 
dence on the link variables. It is not possible to analytically find the minimum of S as 
a function of a single given link, keeping all others (and p{x)) fixed. If / = 0, where the 
scalar degree of freedom is absent, the following algorithm [|ll|, |^ always lowers 5* 



M{U,{x)) 



W^{x) 



\\M{U,{x))-W,{x) 
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We use the same algorithm here and add the prescription for updating the scalar field. 
The definitions of W^{x) and M{Uf^{x)) (specifying the parts of S respectively linear and 
quadratic in f/^(a;)) will be split according to 

M{U,{x)) = M^'\U,{x))+M('\U,{x)) + fM^'\U,{x)), 

W,ix) = + + (x), (19) 

where the index stands for the pure gauge part (k = 0, see ref. the index 1 for 

the K dependent term arising through the modified link equations of motion 0], cmp. 
eqs. (|Tl],|T2|), and the index 2 stands for the part that comes from the scalar equations of 
motion. Using the notation of V^(x) for the 2(n — 1) staples in eq. ([I2|) , we have 

MW(f/^(x)) + M«(f/^(x)) = 2Tr (u,{x)Ul{x)) U,{x) + Gj] Tr (u,{x)V^{xy) V^^{x), 



and 

wj,'\x)=2 y: 



M^'\U,{x))^-^{p\x)+p\x+f,)} Tr (f/,(x)), (20) 



b 



i a^-^i 2r X a X ax 

t JL A* 
W!^'\x) = Kj2p{x+il+a){p{x + a){ ^a]-f\)+p{x+p) pt (| " I )}+p{x)p{x+a) x 

a=^±fj. XX X x + p, X + fi 

(1 ~1) = K^p(x+/i+a)f/„(x){p(x+a)(/-t/^(x + a))[/|(x+/i) + 

X XX ay^ijj, 

p(x + /i) [/^. (a; + a) (^- f^i (a; + /t) ^) } + p(a;) p(x + a) (/- f/^ (x) ) f/^ (x + a) f/] (x + /t) 
Wl,^\x) = Kp{x){'^ J2 p{x+p+a)TT{Ua{x+fi))-p{x+jl){l-2K + 2X{p\x+p)-l))} + 
Kp{x + /i){^ E Tr {Ua{x)) - p(x)(l-2«:+2A(p2(x)-l))}, (21) 

with the unit vectors a, 6 G {±1, ■ ■ ■ , ±n}, and the convention U_a{x) = Ul{x — a). We 
only give the explicit form for Wj^j^^x) and Wj^^x), referring for Wj^^^x) to eq. (19) of 
ref. |]ll|. To implement this algorithm it is useful to point out that Wlf'\x) can be obtained 
by a sum over all links in each staple of U^{x; 0) (see eq. ([T2|) ), with each link Ue replaced 
by the sum over 2Ue{Up — Up), where Up are plaquettes that end at this particular link, not 
overlapping with the original staple. Likewise, Wj^^^x) can be obtained as a sum over all 

links in each staple of f/^(x; 0), with each link Ue{y) replaced by np{y)p{y + £)(/ — Uj{y)). 
Alternatively, one can describe W}^\x)+Wl^\x) by summing over all links in each staple of 
[7^(x;0), replacing each link U,{y) with 2Ut{y){Ul{y)U[{y) - [Ul{y)U',{y)]^}, where tj',{y) 
is defined as f/^d/) in eq. (p!2|), deleting in its sum over staples the one staple that will 
have a link in common with the link (x, x + p). For infinite Higgs self-coupling one puts 
p(x) = 1, and / = to obtain the algorithm of ref. [H]. This is consistent with the fact that 



Wl^^\x)—M^'^\Ufj_{x)) oc K{p{x)dp(x+fi)£+ p{x+ jl)dp(x)£) vanishes when the scalar equations 
of motion are enforced. Note that accidentally ref. M only listed the last of the three terms 



in 



To verify the convergence of this part of the algorithm, we note that 5* changes by the 
following exact amount |jTlI 

6SiU,ix)) = -^Tr {5Ulix){\\MiU,ix)) - W,ix)\\5U,ix) + M{5U,{x))}) . (22) 

For finite values of A, M(f/^(x)) is no longer positive. Nevertheless, for / = 1 and smooth 
configurations (near the continuum limit) one easily sees that M^'^\Up{x)) scales to zero, 
and 6S ^ -112\\6Uf/{g^a^), see ref. Jill (below eq. (24)). 

To complete the description of the algorithm for the general case, we have to specify 
how to update the scalar field. We found that the ordinary cooling, where we replace 
p{x) by ps{p{x)) (eq. ([T5|)) worked well. The apparent reason is that the unstable mode is 
dominated by the gauge part of the energy functional. For large values of this is no 
longer expected to be the case. We have also devised an updating of the scalar field that 
is guaranteed to lower S. Considering only the part Sp^^) that depends on p{x), we find up 
to irrelevant constant factors, 

Spix) = {p{x) [l-2fi:+2A(p2(x)-l)] -p(x)}' + B{x)p{x) + C{x)p^{x), (23) 

where 

B{x) = Kj2{f-^p{x + a) Tr [{I -U^{x))Ul{x; 0)]-[p{x + d){l-2K + 2X{p\x + d)-l)) 

a 

'1 E P{x + M) Tr {U,{x + a))\ Tr (t/„(x))} (24) 



and 



C(^) = ^ Y.{rV{x+h) Tr {I-Ul{x)) + Tr ^(f/,(x))}. (25) 



As before we take a, 6 G {±1, • ■ ■ , and use the convention that U_a{x) = Ul{x—d) and 
f/_a(a;;0) = Ul{x — a;0). Note that S* is a sixth order polynomial in p{x). We will show 
that under very mild conditions 5* is a convex function of p{x). This greatly simplifies the 
problem of minimizing S with respect to p{x), using ordinary Newton- Raphson. Provided 
(aMn)'^ < 12, the second derivative of S with respect to p{x) has a unique minimum at 



Pm{x) = JIps 



[^J^p{x)), with ps defined as in eq. (p!5|). At this minimum 



r _ 16(l-2/.-2A)^-27Ap^(x) + 3A [8(l-2^-2A)p^^(a;)-3p(a:)]^ , 

i^p(-^'Lin- 8(l-2«:-2A) 

As C{x) >0, this is always positive provided 27Ap^(x) < 16(1 —2k — 2A)^, or 

/p(x) \ ' 128(1 - {aMHr/Uf 

[qkvJ 9{aMHy ■ ^ ^ 
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Since p{x)/{6k,v) < p{x)/v, where p{x) is the average over the nearest neighbours, we 
conclude that in all practical cases S is indeed a convex function of p{x). The unique 
minimum of eq. ( p3D is rapidly found by the iteration 

p'ix) = p{x) - sd,^,)S/d%)S, (28) 

where s is a free parameter used to speed up the algorithm (the standard value being 
s = 1). The convexity guarantees that S is always strictly lowered, unless p{x) is already 
at its minimum, like for eq. (|18[). For each sweep one performs both iterations only once 
for each site (one does not gain speed by multiple iterations per site, as the convergence 
of the algorithm is determined by the lowest eigenvalue of the square of the Hessian of the 
energy functional 

Although the algorithm may seem difficult to implement, its main advantage is that it 



is deterministic, with a good understanding of its convergence |TT|. Most importantly, the 
stringent tests that S must always decrease under saddle-point cooling, and the condition 
that for a solution S must vanish to a high degree of accuracy, are guarantees that the 
algorithm was programmed correctly. Also the test for convexity of 5* was never seen 
to be violated after initial ordinary cooling. Testing the algorithm without this initial 
cooling is, even in the absence of the scalar field, not very useful as it tends to get trapped 
in dislocations when starting from a random configuration. This is avoided by ordinary 



cooling due to the choice of positive sign in eq. (^M- 



5 Finite size scaling 

To obtain infinite volume results in the continuum one needs to first extrapolate at a fixed 
volume LMy/ = to the continuum by taking the limit R = 2{aMwY ~^ 0; which 

is achieved by fitting to 

£{MwL, R) = Sspi^iMwL) + Si{MwL)R/2 + S2{MwL)R^/4: + ■■■. (29) 

For small enough lattice spacings this extrapolation can be done accurately. Subsequently 
one extrapolates these continuum results to an infinite volume. The more information 
one has available on the asymptotic behaviour of S{L) the more accurate one can extract 
^s^h = ^sph(oo). Introducing the shifted field ip = (p — (8A)~2M^, we denote by [A, (p) the 
infinite volume solution and by (5lA, Sl^p) the correction due to the periodic boundary 
conditions. The linearized equations of motion are those of non-interacting massive vector 
and scalar fields. For the vector field the linearized equations of motion impose diA'^{x) = 
and the most general rotationally covariant solutions are given by 

Ai{x) = Cw{cos{6)eiajdjK{rMw) + sm{S)M^^eibjeabkdjdkK{rMw)] / y/a^, 

$(x) = CHMHK{rMH)/^/^, K{r) = ^^^^, (30) 

r 

where 6 is non-zero for the deformed sphalerons {Mh > 12Mw) and zero for the ordinary 
sphalerons {Mh < l2Mw)- These functions describe the solution {A,i^) at large distances 
r = ||a;|| —>■ oo. At distances \L > R ^ from the center of the solution, the fields 
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satisfy the linearized equations of motion up to relative errors of the order of e~^^, where 
M is the smallest of the two masses in the problem. In this region the solution can be 
described by periodic copies 

{A{x),ip{x)) = ^ (^(x+nL),$(x+nL)). (31) 

We will now split the energy density V{A + 6lA, (f + Si^p) into V{A, (p) and terms linear 
and quadratic in the shifted fields. Higher order terms are suppressed to (9(e~^*^^/^). To 
this order the term quadratic in the shifted fields, sums with the zeroth order term to 
the energy of the sphaleron in an infinite volume, after integration over the periodic box. 
This is because the dominating contribution for the quadratic term comes from the region 
near the boundary of the periodic box where one can neglect the interactions between the 
copies. To (9(e~^^^/^) all volume dependence is therefore determined by the term linear 
in the shift of the fields, for which we can use the {A, ip) equations of motion, leaving only 
a boundary term 

^,ph(L) = S^^ + d?x d,{5LAt{x)Ff^{x) + 5L^d,^) + Oie-^''"^'^). (32) 

J — L l'2 

The surface integral is evaluated using eq. (0), together with the explicit expressions of 
eq. (^). Each of the six faces of the cube gives the same contribution to the surface 
integral. We extend the integral over one face to the whole plane, at the expense of an 
error of (9(e~^^^^). To this order only the nearest copy will contribute and we can ignore 
the non-linear term in the expression for the field strength. With y = x — iL, one has 

S^^^iL) = S^^ + 6S{L) + 0{e-^''''^), (33) 
6 

aw 

Using dfK{rM) = M'^K{rM), the integrand between curly brackets can be simplified to 

C^cos\5)d.,Kiy)[5adl + ^,^^]K{x) + sm\5)d,K{x)[5,^dl + d.A]K{y) + (34) 

ieakSin{26)[Mwm{y)dkK{x) - d,daK{y)dkdaK{x)/Mw]} + ClM^HK{y)^^K{^^ 

Performing the surface integral one easily sees that the term proportional to sin (25) is a 
total derivative with respect to X2 and X3, whereas at xi = \L the other terms reduce after 
some partial integrations to {K'{r) = dK{r)/dr) 

6S{L) = — f^^^{clMfjK{rMH)K'{rMH) - 2cos{26)C^M^K'{rMw)K{rMw)}. 
ay/ J r ^ J 

(35) 

With = + X2 + xl = \L'^ + and sds = rdr, and the fact that the integrand is a 
total derivative in r, we get the following exact result 

d£{L) = 247rcos(2(^)C^^ -— UnCfj -— . (36) 

aw Mw-I^ aw MhL 

The dimensionless constants 6, Cw and Ch are expected to depend on the ratio Mh/Mw- 
We have thus found the remarkable result that subleading corrections are not powerlike (as 
was assumed in ref. 0), but exponential. With the help of these asymptotic expansions 
we will be able to extract S^^i ^'^ rather high accuracy from our data. 



S£(2xi) = — [ dx2dx3 \Atiy)idiA'^ix)-d^A1ix)) + ^y)di^x)}. 

aw J ^ > 
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6 Results 



As is usual in lattice gauge theories, or for that matter any discretization technique, there 
are two conflicting sources of numerical errors. On the one hand the correlation length 
(1/M) should be much larger than the lattice spacing to minimize lattice artefacts, on the 
other hand it should be much smaller than L = aN to minimize finite size errors. 

For small values of aM^ the electroweak sphaleron tends to develop additional unstable 
modes. There are two reasons due to finite volume effects. The first reason is that the 
rotational invariance will only be approximate such that the energy functional will no longer 
be fiat as a function of the rotational moduli. As saddle-point cooling works irrespective 
of the number of unstable modes, the solution might be attracted to a saddle point with 
additional (usually small) negative eigenvalues of the Hessian. Secondly, the saddle point 
associated to the pure gauge finite volume sphaleron |1T|, obtained by putting k, = 0, 
will be lighter than the electroweak sphaleron for small volumes. At finite values of the 
Higgs self-coupling the pure gauge finite volume sphaleron remains an exact solution by 
putting p{x) = 0. It has an energy 72.605/ (g'^L) + ^M^MfjL^ / g"^. At infinite Higgs self- 
coupling the solution will be deformed (we have no freedom to choose p{x) = to make the 
gauge field massless). In this case the crossing occurs at MwL ^ 2.5. We observed below 
the crossing of these distinct solutions that the electroweak sphaleron acquires additional 
unstable modes. (The other saddle point acquires extra unstable modes for larger volumes. 
Close inspection reveals that the changes do not occur exactly at the crossing.) 

For large values of aMw both translational and rotational invariance will be broken 
by the coarseness of the lattice. This will cause the energy functional to develop spurious 
saddle points and one might get trapped in one with additional negative modes, as for the 
breakdown of rotational invariance due to a finite volume. We typically will choose aM\y 
such that the eigenvalues of the Hessian associated to the approximate zero modes are not 
too big. For finite values of the Higgs self-coupling another feature will cause problems 
at large values of aMw, associated to an enhanced gauge symmetry of the solution. In 
the unitary gauge the energy functional is generally only invariant under global gauge 
rotations. However, suppose that the exact lattice solution will have p(0) = 0, as is true 
in the continuum. It is then easily seen that the hopping term of the energy functional is 
insensitive to all links connected to the origin. The energy functional is therefore invariant 
under a gauge transformation that is non-trivial at x = only, as this does not affect the 
plaquette contribution to the energy. In particular at no expense in energy one can flip 
the sign of the trace of the links connected to the origin. In the way we prepared the 
configurations this will not occur if all links are close to the identity. But at moderately 
large lattice spacing or small volumes p(0) is no longer exactly zero. The hopping term 
now depends weakly on the gauge transformation at the origin. This tends to favour a 
negative value of the trace for only one of the links connected to the origin. (From this 
we also found solutions with the trace of all links positive, with almost identical energies.) 
Initially, a negative value for the trace of one of the links mislead us to believe that we 
were dealing with dislocations. 

Putting all constraints in we found for Mh = oo the window of allowed values to be 
MwL > 2.5, aMw < 0.40, for Mh = Mw the window is MwL > 3.8, aMw < 0.60 and 
for Mh = Ww it is MwL > 4.0, aMw < 0.65. 
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Figure 1 gives the energy density profiles of tlie electroweak spfialeron at eacli of tlie 
tliree Higgs masses. We should not directly use eq. but first average over all directions 
of the links connected to a point x (without affecting the total energy), in order to compute 
the energy density at this point. Note that for = oo the solution is very much more 
peaked in the core region and will have larger lattice artefacts. The behaviour in the 
tail region is similar to the case where M]v = Mh- For Mh = jMw this tail region is 
dominated by the decay of the scalar field. Also plotted in figure 1 is the behaviour of 
p{x)/v for Mh = Mw and Mh = jMy/ at My/L = 4. Because of finite volume effects the 
scalar field does not exactly equal its expectation value at the boundary. Likewise it does 
not quite go to zero at the center, which is also due to finite lattice spacing errors. 




Figure 1: The scalar field (top) and the energy density (bottom) in a plane through the 
center of the electroweak sphalerons for (a) : Mh = oo at MwL = 2.53, (6) : Mh = Mw 
and (c) : Mh = jMw, both at M\yL = 4.0. The energy density is normalized to its peak 
value (respectively 0.093, 0.025 and 0.016M^/aw) and the scalar field p to its expectation 
value V. 

The way we obtained the required configurations was by first constructing a sphaleron 
for the frozen- length Higgs model, starting at = 8. All links at the boundary were 
first put to the identity, which serves the purpose of positioning the solution in the center 
of the lattice and of lifting the energy of the finite volume sphaleron by a considerable 
amount. The latter helps avoid getting trapped in that solution. Centering the energy 
profile will reduce the probability of getting stuck in a saddle point with spurious unstable 
modes due to the breakdown of translational and rotational invariance. We then release 
the frozen boundary condition and compute the Hessian after cooling to verify that we 
have one unstable mode only. This way the maximal energy density occurs at the center 
of a plaquette, see fig. 1. The solutions where the maximum occurs at a lattice point are 
higher in energy. One can now change the lattice spacing in small steps to scan the desired 
range of parameters. For A^ = 12 and 16 the initial configurations were generated from 
the one at A^ = 8, by embedding it in the large lattice (links parallel to the boundary 
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remain constant and those perpendicular to the boundary are put to unity) and adjusting 
the lattice spacing. For = Mh we generated the sphalerons for N = 8 from the 
frozen-length sphaleron (in not too small a volume) by adding the scalar field, set to its 
expectation value v. Varying the lattice spacing in small steps allows one again to scan the 
desired range of parameters. Finally, the sphalerons with Mh = 0.75Mw were generated 
from the ones with M]y = Mh by simply adjusting the parameters. 
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5.09(1) + 13.6(5)^, 
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3.6660(5) 
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3.692 
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-0.23 


4.6000 


3.5490(1) 
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3.6527(4) 


2.308 
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-0.31 
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3.6418(4) 
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-0.11 


3.493 


-0.24 


-0.15 
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-0.24 


-0.17 


3.507 


-0.23 
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4.0000 3.4193(2) 3.4578(2) 3.4703(4) 1.916 

4.4000 3.4078(2) 3.4585(2) 3.4743(3) 1.886 

4.8000 3.3925(1) 3.4584(1) 3.4782(3) 1.934 

5.2000 3.3699(2) 3.4565(2) 3.4807(3) 2.100 



^ = 3.530(3) + 24(4)^ - 12(2)^, ^ = 3.5355 

Table 1: Lattice data for the sphaleron energies S, the negative eigenvalue of the Hessian 
on a 8^ lattice, the fit to the lattice spacing dependence, and volume dependence. We give 
as many digits as we believe to be significant. The variational result is denoted by £y. 



In table 1 we present the results for the sphaleron energies, the negative eigenvalue for 
the N = 8 Hessian, the fit to the lattice spacing dependence (eq. (^)) and to the volume 
dependence (eq. (|36|)). We list the variational results for S^^i &s S^. For the frozen- length 
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Higgs model 0] we have here performed some further coohng down to S < 10~^M^r/(yw 
for = 16 (and one or two orders of magnitude smaller for = 8 and 12), to justify the 
five digit accuracy (estimated errors in the last digit given between brackets). As was to 
be expected, one finds appreciable lattice artefacts for the case Mh = oo. On the other 
hand they are comfortably small for Mh ~ M^. To demonstrate this further we also 
computed at Mw = Mh and N = 8 the energies for aMw = 0.644, 0.663, 0.788 and 0.825 
giving respectively S = 3.493, 3.476, 3.355 and 3.305Mw/(yw- These are solutions with a 
negative trace for one of the links, as described above, which is why we did not use these 
values for the finite size scaling. Nevertheless, it shows that even for these rather coarse 
lattices, the error in the sphaleron energy is only 10%, which was somewhat surprising. 
For these solutions the lattice artefacts are described well by the fit to the lattice spacing 
dependence given in the table (at MwL = 4.8 for Mw = Mh)- 
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3.60 



3.55 
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Figure 2: Continuum extrapolated values for £^sph as a function of the physical volume 
MwL, combined with fits to the finite volume behaviour for (a) : Mh = oo, (b) : Mh = Mw 
and (c) : Mh = Ww- 



Figure 2 compares the fit to eq. (|36|) with the continuum extrapolated lattice data. Our 
results for E^^l very good agreement with the variational analysis [^], particularly 

for Mw = Mh, where we achieve an accuracy of better than .05%. For Mh = oo, a much 
better fit with £^-^ = 5.075{5)Mw/o:w is obtained when dropping the largest-volume data 
point. This was the only case where the energy is lowered significantly as compared to 
ref. 0, seemingly because we are unable to avoid being trapped in saddle points with 
additional unstable modes at A^ > 12. Note that at some point subleading exponential 
corrections will start to become relevant too. For Mh = Mw, dropping the last point gives 
= 3M12{8)Mw/aw, whereas for Mh = \Mw we find 3.b3hMw/aw The values 
of Cfj and cos(2(5)C^ obtained from these fits (cmp. table 1) agree with what one can 
roughly extract from the figures of ref. [Q. 

An alternative method for studying the electroweak sphaleron on the lattice is being 
considered by Ambj0rn and Krasnitz, using the Chern-Simons functional to constrain the 
cooling |12[. It has the advantage of allowing ordinary cooling rather than saddle-point 
cooling and might also be used for computing the energy along the tunnelling path. In the 
continuum the sphaleron has a Chern-Simons number of a half (compared to the trivial 
vacuum). Its disadvantage is that implementing the Chern-Simons functional on a lattice, 
the electroweak sphaleron will only approximately be characterized by a Chern-Simons 
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number of exactly one half. As the tunnelling rate depends exponentially on the sphaleron 
energy, our results might be particularly useful in numerical checks of the semiclassical 
determination [|13| of the tunnelling rate at small temperatures p, Q (for the Abelian 
Higgs model in 1+1 dimensions see ref. 
solution for the sphaleron on a lattice. 



T^), as our method gives the exact saddle-point 
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